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У даній статті розглянуто приклади застосування пакету В для розв'язування рівнянь, які використо- 
вуються для моделювання кібер-фізичних систем. Такі моделі включають диференціальні рівняння, що опи- 
сують неперервну динаміку, а також булеві вирази логіки першого порядку для моделювання дискретних 
подій. Інтегровано диференціальні рівняння за умови, що виконуються умови дискретних подій. У даній 
статті розглянуто приклади використання пакету В дебоїіуе для розв'язання задач у вигляді звичайних дифе- 
ренціальних рівнянь, які розв'язуються лінійними методами. Наведено основні особливості пакету 
К дебоїуе, короткий опис кодів, що дозволяє використовувати його для моделювання різних процесів у нау- 
ці 1 техніці для кібер-фізичних систем. Чисельне моделювання виконувалося в пакеті В для моделі типу 
Лотки-Вольтерри. 
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Вступ 

Стрімкий розвиток науково-техніч- 
ного прогресу, який спостерігаємо в ос- 
танні десятиріччя, потребує появи нових 
методів детекції. Тому в науці, промисло- 
вості та медицині зростає інтерес до кібер- 
фізичних систем (КФС), які є альтернати- 
вою до загальновикористовуваних методів 
вимірювання, що характеризуються пога- 
ною вибірковістю, високою вартістю, по- 
ганою стійкістю, низьким відгуком і пере- 
важно можуть використовуватися лише 
високодосвідченим персоналом. 
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Кібер-фізична система - фізична 
система, яка реалізує інтеграцію обчис- 
лень і фізичних процесів, що відбуваються 
найчастіше у вигляді вбудованих систем і 
мереж для моніторингу та контролю фі- 
зичних процесів у системах зі зворотним 
зв'язком. У таких системах динаміка фі- 
зичних процесів є джерелом інформації 
досліджуваного явища з можливістю кон- 
тролю та розрахунку сигналів керування 
об'єктом |1). 

У КФС дуже часто використовують- 
ся біосенсори, які є новою генерацією сен- 
сорів, що використовують у своїй конст- 
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рукції біологічні матеріали, які надають 
високу вибірковість, селективність, точ- 
ність, дають змогу здійснювати швидкі і 
прості вимірювання |2|. Біосенсори харак- 
теризуються високою ефективністю і ши- 
роко використовуються у харчовій про- 
мисловості |3|, при захисті навколишньо- 
го середовища |4|, в оборонній промисло- 
вості 51, але найчастіше використовують- 
ся у медицині |6-9|, як інструмент для по- 
становки діагнозів. У цілому сімейство 
біосенсорів ділиться на дві частини. Пер- 
ша пов'язана з рівнем рецептора до біоло- 
гічного матеріалу, який використовується 
в його будові. Як рецептори можуть бути 
використані ензим, протеїн, порферін, ан- 
тиген або антитіло. Друга частина біосен- 
сорів обмежена до шару провідника, де 
біологічний ефект перетворюється на ви- 
мірювальний сигнал, який може бути 
електрохімічним, імпедансним, ампермет- 
ричним, оптичним та ін. 

У статті продемонстровано потенці- 
ал використання пакету В для динамічно- 
го моделювання в КФС з метою ефектив- 
ного інтегрування задач, що складаються 
з декількох тисяч змінних стану. 

Постановка проблеми 

КФС належить до високоінтелектуа- 
лізованих інформаційних систем. Вони 
використовують доступний набір інтер- 
фейсів, які дають змогу отримувати швид- 
ку та достовірну інформацію про стан і 
внутрішні дані системи. Згідно з |10| 
КФС, як самоорганізуюча система, потре- 
бує всебічних знань про власну динамічну 
структуру та інфраструктуру загальної 
системи. Для цього необхідно визначити 
типи імуносенсорних пристроїв, врахову- 
ючи їх функціональне застосування. Для 
прикладу, імуносенсори можуть викорис- 
товуватися для оцінки критичних станів 
при серцево-судинних захворюваннях, ве- 
личини інсуліну при вимірюванні величи- 
ни глюкози в крові та виявленні кількіс- 
них показників у деяких фармацевтичих 
сполуках |7-9). 

Саме тому проблему становить роз- 
робка порядку використання пакету 
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В де5оЇуе для математичного моделюван- 
ня динамічних процесів у КФС. 

Аналіз останніх досліджень |і 
публікацій 

Кібер-фізичні системи  ототожню- 
ються з проявом четвертої промислової 
революції, яка відбувається в сучасному 
світі 11). Таким чином, існує також фі- 
зична можливість використання техноло- 
гій «Пиегпеї ої Тбіпєє (Інтернет речей)», 
де необхідно використовувати сигнали від 
давачів і вимірювальних приладів. У 
літературі з'являється все більше публі- 
кацій (12|, які привертають увагу до су- 
часних концепцій та пропонують нові ін- 
новаційні рішення. А. Платцер запропо- 
нував підхід на основі «динамічної логі- 
ки», де описані та проаналізовані кібер- 
фізичні системи |13|, 114). У цих роботах 
використовуються гібридні програми (ГП) 
простою мовою програмування з простою 
семантикою. ГП дають змогу програмісту 
звертатись безпосередньо до дійсних зна- 
чень змінних, які представляють реальні 
величини і визначають їх динаміку. 

У роботі (10) запропоновано загаль- 
ну структуру КФС. При застосуванні цієї 
схеми у випадку імуносенсорів можна ви- 
окремити три напрямки: загальна інфор- 
мація про імуносенсор; вимірювання іму- 
нологічних показників та навиків щодо 
перетворення одиниць та калібрування; 
взаємодія з іншими імуносенсорами. Та- 
ким чином, розглядаються певні методи, 
що дають змогу описати імуносенсор. 

У роботі при дослідженні кібер- 
фізичної імуносенсорної системи вико- 
ристано мову програмування К. Незважа- 
ючи на велике різноманіття мов програ- 
мування, які використовуються при роз- 
робці КФС (Аззетііу, С, СУ, ДЮ, Тама, 
/ауауУстірі, Руфоп, Аа та ін. (151), мова В. 
широко використовується в даний час у 
багатьох галузях, що займаються моделю- 
ванням, машинним навчанням та візуалі- 
зацією даних. 

Мета дослідження 

Запропонувати порядок використан- 
ня пакету ВК де5оЇуе для математичного мо- 
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делювання динамічних процесів у кібер- 
фізичних системах. Навести результати чи- 
сельного моделювання в пакеті ВК де5оіує 
для моделі типу Лотки-Вольтерри. 

Структура імуносенсорів та їх 
характеристики 

Серед великого сімейства біосенсо- 
рів імуносенсори є типовими сенсорами, 
що містять шар рецептора, який чутливий 
і селективний, включаючи імобілізований 
біологічний елемент, наприклад, антиті- 
ло, антиген або хаптен, які є імунологіч- 
ними рецепторами для молекул, які до- 
сліджуються. В імунсорі (імуносенсорі) 
відбувається реакція, яка грунтується на 
взаємодії між антитілом і антигеном або 
маленькими молекулами  (хаптенами). 
Антитіла часто називаються імуноглобу- 
лінами, бо вони є протеїнами, які пов'яза- 
ні з імунною системою. 

Імуноглобуліни використовуються 
імунною системою для ідентифікації та 
нейтралізації чужорідних об'єктів. Вони 
використовують властивості зв'язування 
антигенів. Антигени і антитіла можуть ви- 
користовуватися в шарі рецептора в біо- 
сенсорах. Зменшення властивостей, які 
пов'язані з антитілами під час процесу 
імубілізації антигену, використовуються в 
конструкції шару рецептора, де антитіла 
відіграють функцію аналітів (молекул 
предметного детектування) | 16). 

Молекули, які пов'язані з детекцією, 
забезпечують зв'язування антитіл з антиге- 
нами, утворюючи складні конструкції. При 
цьому між антигенами і антитілами утво- 
рюються дуже сильні зв'язки з константою 
зв'язування Ка - 10" -10"" (171. 

Клітинні біосенсори можуть засто- 
совуватися для кількісної оцінки інфіку- 
вання організму за допомогою певних 
електрохімічних чи оптичних явищ, У ро- 
боті (18) описано клітинний біосенсор, 
який використовує електрохімічну імпе- 
дансну спектроскопію. Даний біосенсор 
призначений для підрахунку людських 
клітин С4-. Область зондування цього 
біосенсора включає в себе електродні пік- 
селі, кожен з яких порівнюється з розмі- 
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ром клітини СІ4-, яка захоплюється пік- 
селями електроду. Вони виявляються 
шляхом спостереження за інформативни- 
ми змінами на пікселі. Стан «Увімкнено» 
або «Вимкнено» електродного пікселя 
вказує на виявлення однієї СІО4-- клітини. 
Таким чином, щоб підрахувати клітини 
СОДА потрібно підсумувати електродні 
пікселі в стані «Увімкнено». 

Цей загальний підхід до кількісного 
виявлення клітин використано для моде- 
лювання імуносенсорної системи, яка за- 
снована на явищі флуоресценції. Імуно- 
сенсори (19) є підгрупою біосенсорів, у 
яких відбувається імунохімічна реакція, 
пов'язана з перетворювачем. Принцип 
роботи усіх імуносенсорів полягає в спе- 
цифічному молекулярному розпізнаванні 
антигенів антитілами для утворення ста- 
більного комплексу. 

Кібер-фізична імуносенсорна сис- 
тема (КФІСС) 

Визначення терміну «Кібер-фізична 
сенсорна система» наведено в (10). Виз- 
начення було введено для промислового 
застосування давачів. Загальне визначен- 
ня КФСС передбачає «більш високий сту- 
пінь поєднання, розподілення системи, 
можливість використовувати / вбудовані 
системи в галузі автоматизації та дотри- 
мання чинних стандартів». Розглянутий 
підхід використано для характеристики 
КФІСС (див. рис. 1), що дає змогу ви- 
конати Його чисельне моделювання. 

Згідно з (10) для визначення КФІСС 
використовуються означення та схеми 
для КФС. КФІСС перетворюють фізично 
виміряні імунологічні показники в циф- 
рову інформацію, яка дає змогу проводи- 
ти обробку сигналів у часі, використо- 
вуючи певні алгоритми. Також є взаємо- 
дія з власними можливостями, вимога- 
ми, внутрішніми даними та внутрішніми 
завданнями з точки зору поширення до 
того ж або вищого рівня ієрархії. 

В основу КФІСС (зовнішній прямо- 
кутник на рис. 1) покладено концепцію 
кібер-фізичної системи (КФС) з врахуван- 
ням особливостей інтелектуальних імуно- 
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сенсорів. З додатковими навиками (пунк- 
тирна лінія на рис. 1) сенсор розширюєть- 
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ше діагностичної інформації про дослід- 
жуваний об'єкт. 


ся до КФІСС. Це дає змогу отримати біль- 
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імуносенсор 


нормований сигнал 


Обробка сигналу 


Електричний 


Імуносенсорна сигнал 
система 


Нормований сигнал 


Електричний 
сигнал 


Масив імунопікселів 
на основі різницевих 


рівнянь 


Вимірювана біологічна величина 


Рис. 1. Функціональна схема КФІСС 


В імуносенсорних пристроях вико- 
ристовуються чотири основні види детек- 
тування: електрохімічний (потенціомет- 
ричний, амперометричний або кондукто- 
метричний (ємнісний), оптичний і термо- 
метричний |19|). Усі типи сенсорів мо- 
жуть використовуватися як прямі (немар- 
ковані) або як непрямі (марковані) імуно- 
сенсори. Прямі сенсори здатні виявляти 
фізичні зміни під час утворення імунного 
комплексу, в той час як непрямі викорис- 
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товують різні рівні генерованого сигналу, 
які дають змогу більш чутливо та універ- 
сально проводити детектування у вимі- 
рювальних системах. 

Найпростіші математичні моделі 
в КФІСС 

На рис.І наведено функціональну 
схему КФІСС для масиву імунопікселів 
на основі різницевих рівнянь, однак у ро- 
боті розглянуто прості математичні моде- 
лі на основі звичайних диференціальних 
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рівнянь (ЗДР) та диференціальних алгеб- 
раїчних рівнянь (ДАРУ. 

ЗДР |20| описують зміну певного 
кількістного показника як функцію часу. 
Вони можуть бути математично представ- 
лені у вигляді: 

' 

УЗ Г/О»Ой), (0) 
де у - диференціальна змінна (характе- 
ризує кількість імунологічного показни- 

/ . . . 
ка), у - похідна, 0 - інша змінна, а Ї - 
незалежна змінна (час). Для розв'язуван- 
ня рівняння (1) задамо початкову умову 
(при І-ї): 


Упр) нс. 

Формальні співвідношення, які на- 
ведені вище, забезпечують явне виражен- 
ня У! як функції у, Х та І. Більш за- 
гальною математичною формою є неяв- 
ний вираз: 

05 СО ум) (2) 

Якщо крім звичайних диференці- 
альних рівнянь на диференціальні змінні 
накладаються деякі алгебраїчні обмежен- 
ня в кожний момент часу: 


0-80.0,1), 0) 
то отримаємо набір диференціальних ал- 
гебраїчних рівнянь (ДАР). Дві попередні 
функції С (формула (2)) та є (формула 
(3)) можна об'єднати до нової функції Б: 

05 ЕС У, 0,1). 4) 

Запропонована функція ЕК буде ви- 
користовуватися в статті. Розв'язування 
ДАР є більш складним, ніж розв'язування 
ЗДР. Наприклад, початкові умови для 
ДАР повинні бути узгодженими. Почат- 
кові значення І, у та у" повинні задо- 


вольняти умові: 


05 КОМИ)» У), о, 1). 5) 

ЗДР та ДАР часто зустрічаються в 

багатьох наукових та інженерних зада- 

чах, наприклад, для моделювання біосен- 

сорних та імуносенсорних інформаційних 

систем |21-22|, електричних кіл, механіч- 
них систем, в рівноважній хімії (231. 
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Більшість ЗДР та ДАР є досить 
складними, тому для спрощення їх аналі- 
тичних орозв'язувань використовуються 
чисельні методи, які дають змогу прово- 
дити обчислення при обмеженій кількості 
значень незалежної змінної 1. 

Спільною властивістю багатьох при- 
кладних програм для чисельних обчис- 
лень є здатність вирішувати «жорсткі» за- 
дачі ЗДР або ДАР |6І. 

До жорстких систем відносяться зав- 
дання хімічної кінетики, нестаціонарні 
процеси в складних системах, що виника- 
ють при вирішенні рівнянь теплопровід- 
ності та дифузії. Жорсткі системи є до- 
сить складними для чисельного розв'язу- 
вання. У таких випадках використовуются 
спеціально розроблені методи, які вимага- 
ють багато часу для обчислень, оскільки 
їм потрібно використовувати дуже малий 
крок дискритизації для задоволення вимог 
до стійкості |71. 

Дуже часто жорсткі системи най- 
більш ефективно вирішуються за допомо- 
гою неявних методів, які вимагають ство- 


рення Якобіана | - | і розв'язування 


У 

системи рівнянь, які пов'язані з цим Яко- 
біаном. За винятком методів Рунге-Кутта, 
всі розв'язування, реалізовані в де5оїує, є 
методами змінного порядку та кроків, які 
використовують формули зворотного ди- 
ференціювання та методи Адамса |З. 

КФІСС включають в себе диферен- 
ціальні рівняння з метою опису неперерв- 
ної динаміки, часові зміни структури 
яких визначаються певними дискретними 
подіями. Тобто, маємо справу із жорстки- 
ми задачами. 

Пакет В де5оїіуе. Характеристика 
пакету В де5оЇуе 

Програмний пакет К використову- 
ється як інструмент для швидкого та 
ефективного моделювання при розробці 
КФІСС. Методи програмної реалізації є 
ефективними, надійними та добре задоку- 
ментованими. Вони включають чотири 
інтегратори з пакету ОРЕРАСК (1. 5ОРЕ, 
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ІЗОРЕЗ, І5ОРА, ТІ5ОРАВ), РУОРДЕ та 
ПАЗРК2.0. Крім того, це набір інтеграто- 
рів Рунге-Кута 1 спеціальних розв'язу- 
вань для ефективного використання ди- 
ференціальних рівнянь у частинних по- 
хідних першого, другого та третього по- 
рядків. Програми розв'язують як жорсткі, 
так і нежорсткі задачі. 

Пакет К де5оіме |9| є наступним 
після пакету оде5оЇуе |10|, який станом 
на сьогодні є частково застарілим. Пакет 
К дебоїме, порівняно з оде8оЇуе, включає 
в себе більш повний набір інтеграторів, 
більш широкий набір параметрів настро- 
ювання і забезпечує більшу продуктив- 
ність. Наприклад, не було можливості 
вказати структуру Якобіанів у звичайно- 
му І5ода (одебоїуе), тоді як така можли- 
вість повністю підтримується дезоіуе. 
Деякі методи інтеграції в дебоїме реалізо- 
вані ефективно, надійно та часто вико- 
ристовуються з відкритим вихідним ко- 
дом. Ці процедури схожі одна на одну і 
добре задокументовані, як у вихідному 
коді, так і в окремих роботах (11). 

Функції І5оде, І5о0де58, І5ода та І5одаг - 
це К-реалізація процедур Когігап з одним 
ім'ям, що належить пакету ОрРЕРАСК 
П2| функції уоде 1 7уоде реалізують 
функції Коггап МОРЕ 1 7МОРЕ (13; 
функція  дазрк реалізує код Когігап 
РА5РК2.0 П14). 

Параметри інтегрування в пакеті 
ВК дебоЇуе 

При виклику процедур інтегрування 
можна вказати багато параметрів з метою 
збереження гнучкості оригінальних кодів. 

Роздільна структура Якобіану має 
вирішальне значення для ефективного ви- 
рішення (помірно) великих систем жорст- 
ких рівнянь. Легкі Якобіани не можуть 
генеруватися набагато швидше, зберігаю- 
чи лише ненульові елементи, що значно 
зменшує вимоги до пам'яті. Одержані 
рівняння можуть бути більш ефективно 
розв'язаними, якщо взяти до уваги скоро- 
чення. Тому користувачі мають можли- 
вість вказати певні властивості Якобіану. 
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За замовчуванням це вважається 
повною матрицею, яка розраховується за 
допомогою розв'язування чисельного 
розмежування (де градієнт функції оці- 
нюється послідовним збуренням змінних 
стану). Щоб скористатися прорідженим 
Якобіаном, можна вказати характерис- 
тику моделі. Таким чином, можна від- 
значити, що Якобіан має стрічкову струк- 
туру (моде, ЧФазрК, Ізоде, І5ода) або вико- 
ристовувати більш загальний нещільний 
Якобіан (І50де8). За винятком методу 
Рунге-Кутта, всі інтегровані методи та- 
кож забезпечують специфікацію аналі- 
тичного Якобіану як варіант, який може 
підвищити продуктивність обчислень. 

Інші важливі компоненти гі! та а 
- відносні та абсолютні допуски, що ви- 
значають помилки контролю. Вони впли- 
вають не тільки на крок, інтегрування, але 
також чисельне розмежування, яке вико- 
ристовується при створенні Якобіанів. 

Короткий опис функцій пакету 
К дебоЇуе 

Усі функції типу І5оде, уоде та дазрК 
використовують диференціювання і змін- 
ним кроком та порядком. Їх доцільно ви- 
користовувати для вирішення жорстких 
завдань (порядку 1-5). Функції І5оде і уоде 
також містять змінні кроки та різні методи 
Адамса (порядку 1-12). 

Функції оде.1Ї), оде.2Р та оде.3Р є 
спеціально розробленими для розв'язу- 
вання рівнянь у частинних похідних, де, 
крім похідної від часу, компоненти також 
змінюються на один, два або три (просто- 
рові) розміри, відповідно. 

Функція І5ода дає змогу автоматич- 
но вибирати жорсткий або нежорсткий 
метод. Це перемикання може відбуватися 
між двома методами під час моделюван- 
ня при зміні жорсткості системи. По за- 
мовчуванні використовується метод оде, 
який найкраще підходить для розв'язу- 
вання простих задач. Функція І5одаг схо- 
жа на Ї50да, але включає в себе метод зна- 
ходження кореня. 

Функції Ізоде та уоде дають змогу 
розв'язувати жорсткі та нежорсткі задачі. 
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Однак користувач має вирішити, який ме- 
тод найбільш підходить для конкретної 
проблеми та обрати відповідний метод 
розв'язування. Код 7уоде є варіантом 
уоде, що розв'язує рівняння за участю 
змінних, які є складними числами. Код 
Ізоде за замовчуванням використовується 
в оде.ІЮ. 

Функція І5оде5 використовується 
для великих жорстких проблем з певними 
взаємодіями між змінними стану, що при- 
зводить до розрідження Якобіанів. При 
використанні І5о0де5 можна досягти знач- 
ної економії часу обчислень. Даний тип 
коду використовується в оде.2Ю і оде.3Р. 

Функція ЧФазрК є єдиним інтеграто- 
ром у пакеті В, який вирішує диференці- 
альні алгебраїчні рівняння індексом нуль 
і один. 

Пакет К також включає розв'язу- 
вання для декількох методів Рунге-Кутта 
(ТК) зі змінним або фіксованим часом. 
Сюди відноситься класичний метод чет- 
вертого порядку Рунге-Кутта та Ейлера 
(ТКА, ешег). 

Крім того, набір коефіцієнтів для 
найбільш поширених методів Рунге-Кутта 
доступний у функції гкМеїроа. (115-171. 

Усі розв'язування повертають ма- 
сив, який містить у своїх колонках зна- 
чення часу (І колонка) і значення всіх 
змінних стану (наступні колонки), за яки- 
ми слідують звичайні вихідні змінні. Цей 
формат особливо підходить для графіч- 
них процедур пакету ВК (наприклад, 
тагріоб). Крім того, додається ріої метод, 
який для моделей з незначною кількістю 
змінних стану дає можливість усі вихідні 
дані вивести на одному графіку. 

Приклад моделювання в пакеті 
К моделі типу Лотки-Вольтерри 

Розглянемо просту модель типу 
Лотки-Вольтерри |18|, що описує взаємо- 
дію «хижак-жертва». Ця модель отримала 
величезну популярність та іноді назива- 
ється класичною моделлю взаємодії по- 
пуляцій хижака і жертви або просто мо- 
деллю «хижак - жертва». Проста матема- 
тична модель добре описує складну біо- 
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логічну, екологічну або економічну сис- 
тему. Прикладом цього служать довго- 
тривалі стосунки між видами хижака 1 
жертви в будь-якій екосистемі. Матема- 
тичні розрахунки зростання популяції ок- 
ремо взятого виду показують, що межі 
щільності популяції можна описати прос- 
тими рівняннями, розв'язання яких дає 
характерну криву чисельності популяції, 
що росте експоненціально, поки вона не- 
велика, а потім вирівнюється, коли вона 
досягає меж можливості екосистеми під- 
тримувати її. Просте продовження цієї 
концепції дозволяє зрозуміти екосистему, 
в якій взаємодіють два види - хижак і 
жертва. 

Взаємодію «хижак-жертва» запише- 
мо у вигляді наступних диференціальних 
рівнянь: 


РЕ) ню 
не , (6) 
райо слі о о 


де Р 1 С - чисельність жертв та хижаків, 
відповідно; г, - коефіцієнт природного 
приросту популяції «жертв» або швид- 
кість зростання чисельності травоїдних за 
відсутності хижаків; /; коефіцієнт 
впливу популяції «хижаків» на чисель- 
ність популяції «жертв» або швидкість, з 
якою зустрічі хижаків з жертвами видаля- 
ють травоїдних з популяції; г, - коефіці- 
єнт природної смертності популяції «хи- 
жаків», або швидкість скорочення чи- 
сельності м'ясоїдних за відсутності траво- 
їдних; К та К,, - характеризують спа- 
дання чисельності жертв і хижаків, відпо- 
відно, внаслідок внутрішньовидової кон- 
курентної боротьби. 

У пакеті В спочатку визначається 
модельова функція, яка називається 
ІУтод0Ю. Необхідно задати вхідні дані: 
поточний час моделювання, значення 
змінних стану та параметри моделі. Ці 
три аргументи завжди повинні бути при- 
сутніми, саме в цьому порядку; інші аргу- 
менти можуть бути додані після них. 
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Во» |Мпод0Ю х- Типспоп(Тіте, 
Зтаге, Раг8) | 

-о о міф(ає П5Кс(951аге, Рагя)), | 
ШтоеіС «-ТІЯР'ЕС 
Стом ПР ««тСЯР'(І-Р/К) 
МогіЄ «тм С 


аР х- Стом ФР - Іеезі С 
аС х- поезі С 7 ДАЕ - МогіС 


гекшгп(15К(С(ЯР, дС))) 
) 
) 

Для параметрів задається ім'я та зна- 
чення (раг5), ініціалізуються змінні стану 
(ушпі) та моменти часу, в які ми хочемо 
вказати їх кількість (йте5). Виходячи з 
цих входів, отримуємо результати моде- 
лювання. Для цього доцільно використати 
функцію інтеграції за замовчуванням оае, 
яка заснована на методі І5ода; кінцевий ре- 
зультат виводиться в плаїгіх, відмітивши 
оці. 


А ЕЕ ЧЕ ЧЕ 


Перша частина матриці (реадй (оио)). 
Матриця ой в першій колонці відліку ча- 


ро 
зесоо реу 
--то ргедаїог 
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су, а в наступних колонках - концентра- 
ції хижаків та жертв. 

К» раг8 х- с(ТІ - 0.2,тС « 1.0,тМ о 
0.2, АЕ «0.5, К - 10) 

К» ушпі « С(Р-1,С -2) 

К» йтез х- 5е4(0, 200, бу - 1) 

К» ргіп(5узіет.те( 

зо оці х- оде(Гипс - ІМпод0Ю, у с 
упші, рагтя - раг5, тез - Ште5))) 

изег зузіет еіарзед 

0.04 0.00 0.04 

ВК» реад(оиії, п - 3) 

йте Р С 

11,7 0 1.000000 2.000000 

2,11 1.626853 1.363283 

13,12 2.438467 1.371156 

В підсумку, модель виводу побудо- 
вано, використовуючи функцію тплабріої, 
пакету В. 

Госка-УМоПетта", Ім - 2) 

К» тагріокоші|, "йте"), оці, 2:31, 
гуре - "1", хіаб - "йте", уіаб - "Сопс", 

- таїп - "роїка-М оПетта", Іма - 2) 

К» Ісєепа("оргієбі", с("ргеу", 
"соп5итег"), сої з І:2, Шу - 1:2) 


6) 


Рис. 2. Результати моделі Лотки-Вольтерри (а) та моделі Лотки-Вольтерри, яка розв'язана 
для стаціонарного стану (ргеу - жертва, ргедаїог - хижак) 


Результати (рис. 2 а, б) показують, 
що після початкових коливань концент- 
рації хижаків і жерт досягають стаціонар- 
ного стану. Це займає 0,04 с (І5ода, дазрК) 
1 0,02 с (П5оде, моде, І5одез) для розв'язу- 
вання цієї моделі. 
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Висновки 

Програмне забезпечення В швидко 
завойовує популярність серед науковців 
та дослідників. Після впровадження паке- 
ту оде5оїме з'явилася можливість вико- 
ристовувати пакет В як інструмент для 


о А.С. Сверстюк 
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розв'язування задач моделювання зви- 
чайними диференціальними. Останні па- 
кети (гооб5оЇме, БурЗ5оїме) пропонують 
для вирішення крайових задач диферен- 
ціальних рівнянь. 

У роботі наведено основні функції 
пакету В. дезоїуе, короткий опис функцій, 
що дає змогу використовувати його для 
моделювання різних процесів у науці та 
техніці при розробці КФІСС. Проведено 
чисельне моделювання в пакеті ВК для мо- 
делі типу Лотки-Вольтерри. 

У подальших дослідження доцільно 
розширити сферу використання програм- 
ного забезпечення В при розробці КФІСС 
на основі диференціальних рівнянь на 
прямокутній та гексагональній решітках. 
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КЕ5ОМЕ 


А.5. 5уег5йШК 

Озіпе В де 5оЇуе расКаге Гог таїе- 
таїйсаї подеїпе ої дупатіс ргосе55е5 іп 
субег-рПузісаї! іптипоз5епзяогу 5у5вет 

Тре аїт ої Фе рарег 15 0 оНег К 50Й- 
хуаге іо 50Їме едиайопя, мбісП аге ц5ед Їог 
тодейпе субег-ріубіса! зузіет5. 5исп то- 
деі5 іпсіцде ФіНегеппа! едиацпопзя дебсгібіпе, 
сопійпиоця дупатіся а5 уме а5 Бооісап ех- 
ргеззіоп5 ої (Де Пг5(-огдег Їобіс Гог 8ітиіа- 
йоп ої дї5сге(е еуепіз. ТПегеїоге ії 15 песе8- 
загу (0 ітіеєртаге ФіНегепиа! едиайопз ргомі- 
дей Фа Аї5сгеіе еуепі5 сопдафійоп5 Бро!4. Пт 
Фі8 агіїсіе уме Фізси85 ехатріе5 ої Ше ц5е ої 
Фе К де5оїуе расКаєе Гог 50Їміпе, ргобіетя 
їп Фе Тогт ої огаїпагу Фійегеппа! едиацопя 
Фаї аге 50Їуед Бу Ппеаг текродз. ТРре таїп 
Геашшге5 ої Ше К де5оіуе расКаєєе, а Бгіеї 
дезсгірйоп ої Ше соде5 аге бїуеп, мбісП 
аПому8 18 їо Бе цеед іо 5ітиіаге уагіоця рго- 
се55е5 їп 5сіепсе апа (есрпоїогу. 

Сопзідегей а 5іпіріе плодеі ої ІГоїКа- 
Мойегга Суре, дез5сгібіпє Ше іпіегасійїоп ої 
"ргедаког-ргеу". А 5ітріе плаїбетайса! то- 
де! де5сгібе5 ммЙ Ше согпріех Біоіовісаї 
зубіетя. Ап ехатріе ої іБіз 15 Ше Іопе-(егт 
теіайоп5рір Бебуееп ргедаког апа ргеу 5ре- 
сіе5 іп апу есо5узіет. Мафетайса! саїси- 
Іапоп5 ої рориіайоп єгоміВЮ ої а 5іпбіе 5ре- 
сіе5 іпдісаїе ас Фе рориіайоп 4еп5іїу П- 
тії5 сап Бе де5сгібед Бу 5ітріе едиайопе, 
хуУро8е 50Гайоп єїмеє а спагасіегі5 с сигуе 
ої рориіайоп 5і7е Шаг ехропепцаПу ягом/5 
ууУрППе 10 15 8плай, апа ШФеп аПпяпей у/реп її 
геаспе5 Ше Птії5 ої Ше есоб5убіет'я абіШу 
іо 5иррогі її. А 8ітріе ехіепзіоп ої (Бі8 соп- 
сері аПом/5 ця іо ппаег5кала Ме есозувіет 
іп уУрбіср бмо 5ресіе5 - ргедаког5 апа ргеу 
іпіегасі. 

Тре рарег ргебепі5 Ше паїп Гипспоп5 
ої Бе Б де5оїЇуе расКаге, а Бгіеї де5сгірйоп 
ої Фе Гипспоп5 Шаг аПомуз її їо Бе ц5ед (0 
5аиіаїе уагіой5 ргосе55е5 іп 5сіепсе апа 
іесппоїогу іп Ше деуеіортепі ої субег- 
рбузіс5 ітпипипозбеп5огу 8у5іетя. Митегісаї 
зитпишайоп у/ає регіогтей їп Фе К расКаєе 
Гог Фе Гоа-Моїетта гуре птодеі. Опе ої 
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Фе адуап(аєєз ої Ше ргоро5ед арргоасі 15 
Фе и5е ої 50ймаге Гог Іасе зузіетзя. Треу 
аге ц5ед їп Ше 5іпиіайоп ої уагіоця субег- 
рбубзісаї! зузіет5 Базед оп ріхе! агпаузє. ТПе 
К де5оіме расКаєе и5е5 єгарріс вгаттаг, 
мрісір сап Бе ицяеад о 4ібріау агеа8, 5исП а5 
ррбазе рогігаїі8, єгарр5 Ша сап гергебепі 
оф сопйпиоц5 апа 4ї5сгеіе срапесе5 іп 
субег-рПузіса! 5у5іет5. 
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